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We review some aspects of the bipolar molecular outflow phenomenon. In particular, 
we compare the morphological properties, energetics and velocity structures of outflows 
from high and low-mass protostars and investigate to what extent a common source model 
can explain outflows from sources of very different luminosities. Many flow properties, 
in particular the CO spatial and velocity structure, are broadly similar across the entire 
luminosity range, although the evidence for jet-entrainment is still less clear cut in massive 
flows than in low-mass systems. We use the correlation of flow momentum deposition rate 
with source luminosity to estimate the ratio / of mass ejection to mass accretion rate. From 
this analysis, it appears that a common driving mechanism could operate across the entire 
luminosity range. However, we stress that for the high-mass YSOs, the detailed physics 
of this mechanism and how the ejected wind/jet entrains ambient material remain to be 
addressed. We also briefly consider the alternative possibility that high-mass outflows can 
be explained by the recently proposed circulation models, and discuss several shortcomings 
of those models. Finally, we survey the current evidence on the nature of the shocks driven 
by YSOs during their pre-main-sequence evolution. 



I. INTRODUCTION 

Stars of all masses undergo energetic, generally bipolar, mass loss during their 
formation. While diagnostics of mass loss using optical and near-infrared line emis- 
sion provide the most dramatic images of the violent birth of stars, molecular flows, 
predominantly traced by their CO emission lines, provide the strongest constraints 
on models of cloud collapse and star formation. This is because molecular flows are 
the massive and predominantly cool reservoir where most of the flow momentum 
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is eventually deposited and so provide a fossil record of the mass loss history of 
the protostar. In contrast, the optical and near-infrared emission arises from hot 
shocked gas which cools in a few years, and hence only traces the currently active 
shocks in the flow. In addition, since there is negligible extinction at millimeter 
wavelengths, even the youngest, deeply embedded objects can be studied. 

Since the last review of this subject in the previous conference proceedings in 
this series (Fukui et al. 1993), many breakthroughs in understanding the nature 
of these objects have been forthcoming, primarily as a result of improvements at 
millimeter and submillimeter wavelength observatories. (We refer the reader in 
particular to volume 1 of the conference series Revista Mexicana de Astronomfa 
y Astroffsica and the proceedings of the IAU Symposium 182). More sensitive 
receivers and focal plane arrays have allowed wide-field imaging of outflows, showing 
that they often extend to many parsecs in size, sometimes even beyond their natal 
cloud boundaries (e.g. Padman et al. 1997). In addition, the availability of mosaiced 
interferometric images of many outflows at 1 — 2" resolution has led to significant 
breakthroughs in understanding the small-scale structure of outflows (e.g. Bachiller 
et al. 1995; Cernicharo and Reipurth 1996; Gueth et al. 1996) and provide strong 
constraints on viable outflow-acceleration models, in particular on the role of jets 
and the nature of the entrainment mechanism (Cabrit et al. 1997; Shu and Shang 
1997). Finally, multi-line studies (Bachiller and Perez Gutierrez 1997) have shown 
that shocks driven by flows chemically process the ISM at a rapid rate, modifying its 
composition and having profound implications for chemical and dynamical models 
of molecular clouds. 

In this chapter, we review some of the recent results on outflows from both low 
and high-mass protostars and focus on a comparison of their physical properties. 
In section II, we summarize the observed flow morphologies and velocity structures, 
and look for evidence that a common mechanism could reproduce the observations 
of all such flows. In section III, the possible entrainment and ejection mechanisms 
for outflows are addressed, and the evidence for a common mechanism is discussed 
based on an analysis of the flow energetics. Finally, in section IV, the evidence for, 
and nature of, shocks in outflows is presented. 

II. OUTFLOW STRUCTURE 

Molecular outflows come in a wide variety of shapes and sizes. This is hardly 
surprising given the likely diversity in the mass and multiplicities of the driving 
protostars, the outflow ages, and molecular environments. Although at least two 
hundred outflows have been cataloged (Wu et al. 1996), there is still no large ho- 
mogeneous sample which has been mapped with good resolution and sensitivity, 
so much of our current understanding is heavily influenced by a few well-studied 
examples. This is even more true for high-mass systems, where perhaps only 10 
or so outflows have been studied in detail, although rapid progress is being made 
using millimeter-wave interferometers to study flows in the often complex environ- 
ment of high-mass star formation (Shepherd et al. 1998). Consequently, our current 
estimates of "typical" outflow properties may be far from accurate. 
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A. Low-Mass Systems 

It is now clear that stellar jets, with speeds of 100-300 km s _1 and densities 
of order 10 3 cm -3 , are responsible for accelerating much of the molecular gas in 
many of the youngest low-mass outflow systems (e.g. Richer et al. 1992; Padman 
and Richer 1994; Bachiller et al. 1995). However, there is also good evidence for 
momentum being deposited into the flows by wind components with wider opening 
angles, and that this component perhaps becomes relatively more powerful as the 
flows age (e.g. Bence et al. 1998). 




Figure 1. The HH211 molecular jet mapped with the Plateau de Bure interferom- 
eter. The main panel shows the high-velocity CO jet in greyscale superposed on 
the lower velocity outflowing gas, which forms a cavity around the jet. The thick 
contours show the 1.3mm continuum emission. The panels top left and lower right 
show the fast and slow CO emission overlaid on the shocked Hi line emission. Data 
taken from McCaughrean and Zinnecker 1996, and Gueth and Guilloteau 1999. 

Some of the first clear evidence for a jet-dominated outflow origin was identified 
in the flows from L1448C (Bachiller et al. 1990; Guilloteau et al. 1992; Dutrey et al. 
1997) and NGC2024-FIR5 (Richer et al. 1989, 1992). In L1448, a spectacular SiO 
jet is seen emanating from the driving source, aligned with the large-scale CO flow 
and unresolved across its width. In the NGC2024-FIR5 flow, the fastest gas is seen 
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to form an elongated jet-like feature on the axis of nested shells of lower velocity 
gas. The collimation ratio q, defined as the width of the flow to the distance to the 
driving source, is as high as 30 for the high velocity (30 km s _1 ) gas in this source, 
but only 4 or so for the lower velocity (5 km s _1 ) envelope. Many other flows are now 
known which show this structure, with high-velocity elongated components tracing 
jet-like activity lying inside cavities of lower velocity gas; examples include L1157 
(Gueth et al. 1996), HH111 (Cernicharo and Reipurth 1996; Nagar et al. 1997) and 
HH211 (Gueth and Guilloteau 1999). All of these flows appear to be driven by very 
young, low-mass objects, based on their low luminosities, and non-detection at even 
infrared wavelengths. The CO flow from HH211 (Fig. 1) is perhaps the most striking 
image to date of this phenomenon, showing an unresolved CO jet with high velocity 
gas (v > 10 km s _1 with no correction for inclination) lying within an ovoid cavity 
of slower moving gas (< 10km s _1 ). HH211 is probably also one of the youngest 
low-mass outflows known, having a dynamical age ra = 0.07pc/10km s _1 = 7000 
years; if the source is inclined close to the plane of the sky, as is expected given the 
clear separation of red and blue outflow lobes and the relatively modest projected 
flow speeds, the true dynamical age could well be a factor of 5 or so lower. It appears 
that HH211 is at the very start of its main accretion phase, and this perhaps explains 
the relative simplicity of the outflow structure. Note that the full opening angle at 
the base of the flow is only 22°. 

It is natural to identify these so-called "molecular jets" as the deeply-embedded 
counterparts of the Herbig-Haro (HH) jets seen in less obscured systems and as the 
neutral counterparts of the ionized jets seen in the radio (see chapters by Eisloffel 
et al., and by Hartigan et al., this volume). Important work by Raga (1991) and 
Hartigan et al. (1994) demonstrated that HH jets were denser and hence more 
powerful than initial estimates suggested (e.g. Mundt et al. 1987), so that the 
total momentum flux in HH jets integrated over the lifetime of a typical source is 
sufficient to drive the observed CO flows (e.g. Richer et al. 1992; Mitchell et al. 
1994). This unified picture of jet-driven outflows is entirely consistent with the 
observed CO structures discussed above. (Masson and Chernin 1993; Cabrit et al. 
1997; Smith et al. 1997; Gueth and Guilloteau 1999). However, we caution that 
the actual composition of the driving jet, whether primarily atomic or molecular, is 
unknown. It is still unclear whether the CO jets seen in sources such as HH211 and 
NGC2024-FIR5 arises (1) from the body of a jet where molecules have formed in 
the gas phase (Glassgold et al. 1989; Smith et al. 1997); (2) from molecules formed 
in the post-shock region of shocks in the jet; or (3) from ambient molecular gas 
turbulently entrained along the jet's edge and at internal working surfaces (Raga 
et al. 1993; Taylor and Raga 1995). 
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Figure 2. CO 2-1 image (in greyscale) of the RNO 43 molecular outflow (Bence at 
al. 1996); the driving FIR source is at the origin, marked with a cross. The solid 
white patches show the extent of the Ha emission. Note the large flow extent and 
correlation of the Ha emission with the CO hotspots. The panel to the right shows 
a detail from the northernmost Ha emission patch, showing the bow-shock shaped 
Ha structure in greyscale, and CO 2-1 contours overlaid. 

RN043-FIR is another low-luminosity outflow source, and is probably a flow in 
middle-age (Bence et al. 1996; Cabrit et al. 1988): the driving source, with L* = 6 
L©, is a heavily embedded Class protostar invisible even at 2/im, but the flow 
extends over a total size of 3 pc. The outflow axis is close to the plane of the sky, 
so the image of CO integrated intensity shown in Fig. 2 shows both blue and red- 
shifted sides of the flow. The dynamical age is ra ~ 1.5pc/10km s _1 = 1.5 x 10 5 
years which is an upper limit due to the small but unknown inclination angle. The 
CO flow is much more complex than HH211 and this most likely reflects the clumpy 
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nature of the molecular gas through which the driving jet has passed; nonetheless, 
a very approximate S-shape symmetry can be seen about the driving source, both 
in the spatial and velocity data, suggesting that the jet direction has changed over 
time (Bence et al. 1996). A precessing jet model, with a full cone angle of 27° and 
a precession timescale of order 3 x 10 4 years, provides a reasonable description of 
the overall flow properties (Bence et al. 1996); such precession could be driven by 
a binary companion in an inclined orbit. 

The evidence that this flow is jet-driven even at a distance of 1.5 pc from the 
driving source is demonstrated by the Ha images of the object, which show strong 
emission coincident with the brightest CO points (see Fig. 2). In addition, at the 
northernmost CO feature in the flow, bright Ha coincident with a bow-shock shaped 
CO feature suggests that the whole of the RN043 outflow can be explained by the 
gradual sweeping up of a clumpy molecular cloud by a powerful stellar jet whose 
ejection axis varies slowly with time. It is also interesting to note that parsec-scale 
flows such as RN043 are the natural counterparts to the parsec-scale Herbig-Haro 
flows (Reipurth et al. 1997) seen in wide-field optical imaging: if there is molecular 
gas in the jet's path, then flows such as RN043 result, whereas jets in essentially 
empty space such as HH34 (Devine et al. 1997) show only optical emission. This 
points strongly to the jets being primarily atomic in composition. However, in the 
HH111 outflow system, CO "bullets" associated with the optical jet and having 
similar velocities to the optical gas have been detected far beyond the molecular 
cloud boundary (Cernicharo and Reipurth 1996). This suggests that in some cases 
the jets may have a molecular component. 

Not all low-mass flows are jet-dominated: many of the older flows show CO 
emission dominated by low-velocity cavities with little evidence for elongated high- 
velocity CO jet features. Examples include L43 (Bence et al. 1998), L1551 (Moriarty- 
Schieven et al. 1988) and B5 (Velusamy and Langer 1998). These flows are typically 
10 5-6 years old, and associated with nebulosities visible in the optical or at 2 ^m, 
suggesting they are Class I or II protostars, older than the Class objects responsi- 
ble for the HH211, RN043 and NGC2024-FIR5 outflows. The lack of high-velocity 
CO and obvious jet-like features in these flows, in contrast to their younger coun- 
terparts such as HH211, and the presence of much wider cavities at the base of the 
flows, suggests that the jet power has declined over time, and that a wider opening 
angle wind is now primarily responsible for driving the outflow. However, even in 
these older sources there is usually some evidence for weak jet activity: L1551 has 
an optical jet and extended Herbig-Haro emission apparently on one of its cavity 
walls (Davis et al. 1995; Fridlund and Liseau 1998), as well as fast CO emission sug- 
gestive of a jet-origin (Bachiller et al. 1994), and B5 has optical jets stretching 2.2 
pc away from the driving source. However, the L43 flow shows no signs of shocks, 
jets or very fast CO, and this may be a true "coasting" flow, which is no longer 
being accelerated by a stellar wind or jet. 

The recent interferometric images of the B5 outflow (Velusamy and Langer 
1998) reveal a beautiful example of wide, hollow cavities at the base of the out- 
flow, much like the reflection nebulosities seen in the near-infrared in many of these 
systems; in order to reproduce the clearly separated red- and blue-shifted emission 
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lobes, the CO must be flowing along these cavity walls, presumably being acceler- 
ated by a poorly-collimated radial wind from the star. The entire outflow driven 
by B5 has a dynamical age of 10 6 years, and a collimation factor of about 5. The 
cavity opening angle at the protostar is extremely large, in the range 90 — 125°, 
and Velusamy and Langer (1998) suggest that if this angle further broadens with 
time it may ultimately cut off the accretion flow. This idea is consistent with most 
of the available data: the young flows such as HH211, RN043 and L1157 have 
opening angles less than 30°or so, while the older flows such as L1551, L43 and B5 
are significantly greater than 90°. With maps at good resolution of a larger sample 
of outflows, it will be possible to test the hypothesis that flow opening angle is a 
measure of the source age. 

B. High-Mass Systems 

Our understanding of massive flows is beginning to change because we are 
starting to find a few isolated systems that can be studied in depth. However, we 
are still observationally biased toward older flows that are easier to identify and 
study. This bias is likely to diminish in the future as more massive young flows are 
identified (e.g. Cesaroni et al. 1997; Molinari et al. 1998; Zhang et al. 1998). 

Most luminous YSOs have relatively wide-opening angle outflows as defined 
by their CO morphology. Although the statistics are poor because few massive 
flows have been studied with sufficient resolution to adequately determine the mor- 
phology, collimation factors q for 7 well-mapped flows produced by YSOs with 
Lboi > 10 3 L range from 1 to 1.8 (NGC 7538 IRS1: Kameya et al. 1989; HH 80-81: 
Yamashitaet al. 1989; NGC 7538 IRS9: Mitchell et al. 1991; GL 490: Mitchell et al. 
1995; Ori A : Chernin and Wright 1996; W75N: Davis et al. 1998; G192.16: Shep- 
herd et al. 1998). The dynamical time scales for these outflows range from 750 years 
to ~ 2 x 10 5 years and there is no obvious dependence of flow collimation on age. In 
comparison, collimation factors in low-mass outflows range from ~ 1 to as high as 
10 with a typical value being ~ 2 or 3 (Fukui et al. 1993 and references therein). It 
appears that more luminous YSOs do not in general produce very well-collimated 
CO outflows and this result is independent of outflow age, unlike outflows produced 
by low-luminosity YSOs. This may be due to the fact that outflows from lumi- 
nous YSOs tend to break free of their molecular cloud core at a very early stage in 
the outflow process. Hence, massive molecular flows are frequently the truncated 
base of a much larger outflow that extends well beyond the cloud boundaries (e.g. 
HH 80-81: Yamashita et al. 1989; DR21: Russell et al. 1992; G192.16: Devine et 
al. 1999). 

The outflow from Orion A is perhaps the best known and best studied high-mass 
outflow system. The driving source is believed to be an O star, and the estimated 
age of the flow is ~ 750 years which makes this one of the youngest known outflows. 
The flow differs significantly from those produced by low-luminosity sources and 
represents a spectacular example of the outflow phenomenon. The CO outflow is 
poorly collimated at all velocities and spatial resolutions and the H2 emission is 
dispersed into a broad fan shape that is unlike any other known outflow. Its mor- 
phology is highly suggestive of an almost isotropic, explosive origin. McCaughrean 
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and Mac Low (1997) model the H 2 bullets as a fragmented stellar wind bubble 
using the fragmentation model of Stone et al. (1995) and suggest that the bullets 
are caused by several young sources within the BN-KL cluster. Chernin and Wright 
(1996) argue that the flow is driven by a single massive YSO, source I. The esti- 
mated opening angle, corrected for inclination, is approximately 60° for the blue 
lobe and 120° for the red lobe. 




Figure 3. The G192.16 outflow mapped with the OVRO interferometer. Contours 
of redshifted CO (thick lines) and blueshifted CO (thin lines) delineate the bipolar 
flow emanating from a dense core traced by the central dust continuum emission 
represented in greyscale. 

The outflow from G192.16 is perhaps more typical of outflows from B stars. 
Figure 3 shows an interferometric CO J=l-0 image of the outflow together with 
a closeup of the 3 mm continuum emission showing a flattened distribution of hot 
dust (Shepherd et al. 1998; Shepherd and Kurtz 1999). The 3 pc long molecular 
flow represents the truncated base of a much larger flow identified by Ha and [SII] 
emission that extends almost 5 pc from the YSO (Devine et al. 1998). The CO flow 
is ~ 10 5 years old and appears to be driven by an early B star. Despite the very 
different masses and energies involved, the CO morphology looks very similar to 
L1551-IRS5: the extent of the G192.16 CO outflow is 2.6 x 0.8 pc, approximately 
twice that of L1551, but the mass in the outflow 95 M ) is much greater. The 
opening angle at the base of the flow is ~ 90°. 

The present lack of well-collimated CO outflows with q > 3 from YSOs with 
Lboi > 10 3 L© does not mean that jets and well-collimated structures are not present 
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in these massive sources. For example, the central source in HH 80-81 (Lboi ~ 2xl0 4 
Lq) powers the largest known Herbig-Haro jet with a total projected length of 5.3 
pc, assuming a distance of 1.7 kpc (Marti et al. 1993, 1995 and references therein). 
However, the CO flow appears poorly collimated with q ~ 1 when mapped at 
moderate resolution (Yamashita et al. 1989). Also, the biconical thermal radio 
jet from Cepheus A HW2 (Lboi ~ 10 4 L©) appears to be responsible for at least 
part of the complicated molecular flow seen in CO and shock-enhanced species such 
as H2, SiO, and SO (e.g., Doyon and Nadeau 1988, Martin-Pintado et al. 1992, 
Hughes 1993, Torrelles et al. 1993, Rodriguez et al. 1994, Rodriguez 1995, Garay 
et al. 1996, Hartigan et al. 1996, Narayanan and Walker 1996). The HH 80-81 
and Cepheus A HW2 systems demonstrate that high-luminosity YSOs can produce 
well-collimated jets like those found in association with less luminous stars, even 
though the CO flow may appear chaotic or poorly collimated. Other examples of 
possible jets in massive outflows include IRAS 20126 and W75N IRS1. IRAS 20126 
(L bo i ~ 1.3 x 1O 4 L ) appears to drive a compact jet seen in SiO and H 2 (Cesaroni 
et al. 1997). W75N IRS1 (L bo i ~ 1.4 x 1O 5 L ) shows H 2 2.12^/m shock-excited 
emission at the end and sides of the CO lobe, with a morphology and emission 
characteristics highly suggestive of a jet bowshock (Davis et al. 1998). However, 
this "bowshock" is 0.3 pc wide, i.e. 30 times larger than the H 2 2.12(Um bowshock at 
the end of the HH211 flow (cf. Fig. 1). It is not fully clear whether such wide bows 
are created by protostellar jets, or by a low collimation wind component. Scaling up 
from current hydrodynamical jet simulations (e.g. Suttner et al. 1997), one would 
need a jet radius of ~ 0.03 pc at 1.3 pc from the star, hence a jet opening angle ~ 
2.6°. 



C. Outflow Velocity Structure 

From the above discussion, we conclude that jet activity, bow shocks, and CO 
cavities are common to outflows from low and high-mass systems. There is some 
evidence that high-mass systems are less well collimated than low-mass ones, but 
this may be due to selection effects — young, high-mass systems with small opening 
angles may simply be missing from the small sample currently known (although the 
Orion outflow does appear to be very young) . This conclusion suggests it is sensible 
to consider the possibility that a common driving mechanism is responsible for all 
outflows. 

Several authors have noted that molecular flows seem to be characterized by a 
power-law dependence of flow mass Mco(v ) as a function of velocity. The power-law 
exponent 7 (where Mco(v) oc v 1 ) is typically ~ —1.8 for most low-mass outflows, 
although the slope often steepens at velocities greater than 10 km s _1 from vlsr 
(e.g., Masson and Chernin 1992; Rodriguez et al. 1982; Stahler 1994; Chandler et 
al. 1996; Lada and Fich 1996; Gibb, personal communication). 
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Figure 4- The slope 7 of the mass spectrum M(v) plotted as a function of (a) source 
bolometric luminosity and (b) flow dynamical age t^. Triangles represent 7 for gas 
with projected speeds less than 10 km s -1 relative to the source, and squares are 
for gas moving at more than 10 km The dashed, horizontal line separates the 
low and high-velocity 7 's for clarity. Sources with Lboi < 10 3 Tq are plotted with 
open symbols, those with Lboi > 10 3 L with filled symbols. Representative error 
bars are displayed in the lower left corner of each plot. The solid line in Fig. 4b is 
a linear least squares fit (slope —1.7 ± 0.6 ^ to high-velocity 7's in luminous sources 
versus log(t dy n)- The sources plotted here are VLA1623, IRAS03282, LI448-C, 
L1551-IRS5, NGC2071-IRS1, and Ori A IRC2 (Cabrit and Bertout 1992); TMC-1 
and TMC-1A (Chandler et al. 1996); L379-IRS1-S (Kelly and MacDonald 1996); 
NGC2264G (Lada and Fich 1996); G5. 89-0. 39 (Acord et al. 1991). CephE (Smith 
et al. 1997) HH251-254, NGC7538-IRS9, W75N-IRS1 and NGC7538-IRS1 (Davis 
et al. 1998); G192.16 (Shepherd et al. 1998); and HH26IR, LBS17-H, G35. 2-0.74, 
HH25MMS (Gibb, personal communication). 

Figure 4 plots (a) 7 versus Lboi an d (b) 7 versus the dynamical time scale 
tdyn for a new compilation of 22 sources with luminosities ranging from 0.58 L to 
3 x 10 5 L Q . Triangles represent slopes derived from gas moving less than 10 km s _1 
relative to vlsr while squares represent slopes derived from gas moving more than 
10 km s _1 relative to vlsr- Sources with Lboi < 10 3 L are plotted with open 
symbols while those with Lboi > 10 3 L are plotted with filled symbols. Both 
well-collimated and poorly-collimated outflows are represented in the sample. 

The most striking result from Fig. 4a is that 7's for low-velocity gas are similar in 
sources of all luminosities. This suggests that a common gas acceleration mechanism 
may operate over nearly six decades in Lboi- In addition, there is a clear separation 
between 7's in high and low-velocity gas which supports the interpretation that 
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there are often two distinct outflow velocity components, perhaps corresponding to 
a recently accelerated component and a slower, coasting component. 

Hydrodynamic simulations of jet-driven outflows from low-luminosity YSOs 
predict such a change of slope at high velocity. They also predict that 7 should 
steepen over time possibly due to the collection of a reservoir of low-velocity gas 
(Smith et al. 1997). Figure 4b reveals marginal evidence that the mass spectrum in 
flows from luminous YSOs does becomes steeper with time, both in the low and high 
velocity ranges. The solid line in Fig. 4b is a linear least squares fit (slope — 1.7±0.6) 
to high- velocity 7's in luminous sources versus log(td yn )- There is no indication of 
time evolution of the mass spectrum slopes in outflows from low-luminosity sources. 

The decrease of 7 with time in more luminous sources may be due to a dif- 
ference in the driving mechanism, or it may simply be more prominent because 
the mass outflow rate is several orders of magnitude greater than in outflows from 
low-luminosity YSOs (thus allowing more precise determination of 7) and because 
their flow ages cover a broader range, from 750 to 2 x 10 5 yrs. 

III. TESTS OF PROTOSTELLAR WIND AND ACCRETION MOD- 
ELS 

A. Flow Energetics 

It is well known that outflow energetics correlate reasonably well with L ho \ over 
the entire observed luminosity range. In Fig. 5a, we show a recent compilation of 
the mean momentum deposition rate Fco as a function of bolometric luminosity 
of the driving star. It must be remembered that -Fco is the time-averaged force 
required to drive the CO outflow: Fco = Mqo vqo/t&, where we assume ra is a 
good approximation to the flow age. If the CO-emitting material is accelerated by 
a separate stellar wind (or jet), the wind momentum flux F w may be quite different 
from Fco 5 depending on the nature of the wind-cloud interaction. There are clearly 
large uncertainties in the measured flow properties, primarily due to difficulties in 
estimating the inclination angle of the outflow, the optical depth and excitation of 
the CO, and the correctness of the assumption that the dynamical timescale is a 
good estimate of the flow age (Cabrit and Bertout 1992; Padman et al. 1997). How- 
ever, as seen in Fig. 5a, the correlations are roughly consistent with a single power 
law (here of slope ~ 0.7) across the full range of source luminosities. It has been 
suggested that this correlation argues for a common entrainment and/or driving 
mechanism for molecular flows of all masses; but this is by no means a compelling 
argument, given that we would surely expect most physically reasonable outflow 
mechanisms to generate more powerful winds if the source mass and luminosity are 
increased. 

Regardless of the details of how the stellar wind or jet entrains material, if the 
shock cooling times are short compared to flow dynamical timescales (as we expect 
for wind speeds less than 300 km s : Dyson [1984], or if there is efficient mixing at 
the wind/molecular gas interface: Shu et al. [1991]), then the wind and molecular 
flow momenta will be equal. This is often called the momentum-conserving limit. 
Then, molecular outflows represent a good opportunity to test proposed protostellar 
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ejection mechanisms. In particular we show in this section that they can be used 
to estimate the ratio / of ejection to accretion rates that would be necessary if the 
wind is accretion-powered. 

We first recall that optically thin, line-driven radiative winds (such as those 
present in main-sequence O and B stars) are insufficient to drive the flows if the 
flows are momentum driven. The solid line in Fig. 5a shows the maximum momen- 
tum flux available in stellar photons in the single-scattering limit, L\, \jc. It falls 
short of the observed amount in molecular flows by one to three orders of magnitude. 
If flow lifetimes have been underestimated, the discrepancy is reduced, but not by 
a sufficient amount. In principle, higher momentum flux rates could be reached 
with multiple scattering. Wolf-Rayet stars (Lboi ~ a ^ ew 10 5 L©) have a wind 
force reaching 20 — 50Lb o i/c, similar to the momentum rate in molecular flows from 
sources of comparable luminosity (cf. Fig. 5a). Such high values are attributed to 
multiple scattering of each photon by many lines closely spaced in frequency (e.g. 
Gayley et al. 1995 and refs. therein). However, excitation conditions in protostellar 
winds are very different from those in hot, ionized Wolf-Rayet winds. The dashed 
line in Fig. 5a shows the typical momentum flux in the ionized component of proto- 
stellar winds, inferred from recombination lines or radio continuum data (Panagia 
1991): the values lie a factor of 10 below Fco, implying that the driving winds 
must be 90% neutral. If dust grains instead provide the dominant opacity source in 
protostellar winds, comparison with winds from cool giants and supergiants might 
be more appropriate. Their wind velocities are ~ 5-30 km s _1 and mass- loss rates 
do not exceed 10~ 4 M yr _1 , hence F w < 0.5 — 3Lboi/c for a typical Lboi ~5x 10 4 
L . Calculations by Netzer and Elitzur (1993) show that 10 times larger mass- loss 
rates could in principle be achieved in oxygen-rich stars, where silicates dominate 
the opacity curve. The wind force could then become comparable to the molecu- 
lar outflow momentum rate for Lboi ~ 5 x 10 4 Lq. Hence it is just possible that 
luminous stars (L* > 5 x 10 4 L©) could drive their flows by radiative acceleration 
if dust opacity plays a significant role; further work is needed to investigate if such 
models are viable. In lower luminosity sources, however, the opacity required to lift 
material above escape speeds largely exceeds typical values for circumstellar dust. 
Therefore, molecular outflow sources with L* < 5 x 10 4 Lq must possess an efficient 
non-radiative outward momentum source. 

It thus appears most likely that both low and high-mass systems possess an 
efficient non-radiative wind-generation mechanism in their embedded protostellar 
phase. If energetic bipolar winds are the chief means of angular momentum loss 
during the main accretion phase for stars of all mass, this is not surprising. However, 
the details of the wind ejection mechanism could differ; in particular massive ac- 
creting stars are likely to have thinner convective layers and probably rotate faster, 
so that the magnetic field and accretion geometry close to the star may be very 
different from that in low-mass stars. 

There exist quite a number of efficient accretion-powered wind mechanisms 
from the stellar surface, disk, or disk-magnetosphere boundary. The most efficient 
models use a strong magnetic field in the star or disk to drive the wind and to 
carry off angular momentum from the accreting gas (see e.g. reviews by Konigl and 
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Pudritz, and Shu et al., this volume). This wind is further collimated by magnetic 
or hydrodynamic processes (e.g. Mellema et al. 1997), generating a high-Mach 
number wind and/or jet with a speed of order 200-800 km s _1 . A fraction / < 1 of 
the accretion flow M a is ejected in the wind: M w = f M a . 

A different scenario recently explored by Fiege and Henriksen (1996a, 1996b) 
is that molecular flows are not predominantly swept-up by an underlying wind, 
but represent infalling gas that has been deflected into polar streams by magnetic 
forces. Only a small fraction of the infalling gas actually reaches the star to produce 
accretion luminosity. In that case, / > 1. This model has been invoked in particular 
to explain the large flow masses Mqo ~ 10M* observed in flows from luminous 
sources. In the following we use molecular flow observations to set constraints on 
these two classes of proposed models. 

B. The Ejection/ Accretion Ratio in Protostellar Winds 

Two simple, independent methods have been used in the literature to estimate 
/, assuming that the driving mechanism is steady over the source lifetime. First, in 
low-mass protostars where the luminosity is accretion-dominated, it is possible to 
use the observed correlation of flow force with Lboi (Fig. 5a) to derive the ejection 
fraction /. The source bolometric luminosity is Lboi = GM*M a /R* = M a VK 2 , 
where vk is the Keplerian speed at the stellar surface. The flow and wind force (as- 
suming momentum conservation) is Fqo = f M a v w . Hence Lco/Lboi = / v w /vk 2 - 
A value / ~ 0.1 is inferred in both Class and Class I low-luminosity objects 
(Bontemps et al. 1996). Alternatively, if one estimates M* from Lboi via the ZAMS 
relationship (which is probably only valid for sources with Lboi> 10 3 L ), one can 
use the accumulated flow momentum to infer /, using: Leo = v w M w = f M*v w . 
Values of / ranging from 0.1-1 are inferred for a wind speed of 150 km s -1 (Masson 
and Chernin 1994; Shepherd et al. 1996). However, we point out that wind speeds 
are unlikely to remain constant over the whole Lboi range. There is evidence for 
higher wind velocities in luminous flow sources: for example, proper motions up to 
1400 km s _1 are seen in HH80-81 (Marti et al. 1995) and Z CMa shows optical jet 
emission with speeds up to 650 km s -1 (Poetzel et al. 1989). These are significantly 
higher than the 100-200 km s _1 typically seen in low-mass sources. 

To re-examine this issue in a homogeneous way over the whole luminosity range, 
we plot in Fig. 5b the values of fv w /vx obtained using a new combination of the 
above two methods (Cabrit and Shepherd, in preparation). The solid circles are 
derived assuming the luminosity is accretion dominated, while the open circles (for 
sources more luminous than 1O 3 L0) assume the ZAMS relationship given above. 
Typically, Vk ranges from 100 km s _1 in low-luminosity sources to 800 km s _1 in 
high-luminosity sources. A rather constant value fv w /vK ~ 0.3 seems to hold over 
the whole range of Lboi- This value is in line with both popular MHD ejection 
models. In the X-wind model (Shu et al. 1994, and this volume), the wind is 
launched close to the stellar surface (y w ~ %) and a large fraction of the accreting 
gas is ejected (/ ~ 0.3). In the self-similar disk- wind models (Ferreira 1997; Konigl 
and Pudritz, this volume), less material is ejected (/ ~ 0.03) but the long magnetic 
lever arm accelerates it to many times the Keplerian speed (v w ~ 10vk)- Thus 
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we conclude that the energetics of the flows over the entire luminosity range are 
broadly consistent with a unified MHD ejection model for all flow luminosities, but 
we reiterate that given the very different physics involved around high and low-mass 
protostars, the details of such a model for high-mass systems are still unclear. 




log Lbol (Lo) log Lbol (Lo) 



Figure 5. (a) The average momentum flux in the CO flows as a function of source 
luminosity. The solid line shows the force available in stellar photons (assuming 
single scattering), and the dashed line the force available from the ionized wind 
components (Panagia 1991). (b) The factor fv w /vK derived for the same set of 
sources is plotted as solid circles (small open circles assume that luminous sources 
are on the ZAMS). Large open circles show the required (U(9)) factors for the same 
objects if the circulation models are applicable (see text for details). 

We stress again that this plot assumes perfect momentum conservation in the 
wind/flow interaction. There could be strong deviations from this key assumption. 
First, we should keep in mind the possibility that massive flows enter the energy- 
driven regime: for high wind speeds, the gas cooling behind the shock will be slow, 
and the snowplow or momentum-conserving flow will turn into an energy driven one. 
In this case, the shell momentum can exceed the momentum in the wind itself by a 
factor of order v w /vco (Cabrit and Bertout 1992; Dyson 1994) so reducing markedly 
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the momentum requirement of the driving wind. Of course, there are objections to 
energy driven flows, in particular their inability to reproduce the bipolar velocity 
fields of many flows (Masson and Chernin 1992); but given the apparently poorer 
collimation of high-mass outflows this issue should perhaps be re-examined. Second, 
if the flow is entrained in a jet bow-shock, the efficiency of momentum transfer will 
depend on the ratio of ambient to jet density: it will only be close to 1 if the jet 
is less dense than the ambient medium. Highly overdense jets will pierce through 
the cloud without depositing much of their momentum (Chernin et al. 1994), and 
in that case the ratio fv w /vx plotted in Fig. 5b would have to be increased. We 
conclude that for current protostellar jet models to apply, we have the additional 
condition that most of the jet momentum must be in a component that is not 
significantly denser than the ambient molecular cloud on scales of 0.1-lpc. 

The above analysis also provides important constraints on accretion rates in 
protostellar objects of various masses: if fv w /vK ~ 0.3, then the values of Fco m 
Fig. 5a show that M a = 3Fqo/vk must range from a few 10~ 6 M yr _1 at L bo i~ 
ILq to a few 10~ 3 M yr _1 at Lt>oi~ 10 5 L©. Then, protostellar sources would not 
be characterized by a single infall rate across the whole stellar mass range, and 
massive stars would form with much higher infall rates than low-mass stars (Cabrit 
and Shepherd, in preparation). 

The wealth of data available on low-mass systems also allows one to break the 
samples down by estimated age, and so look for evolution of outflow properties with 
stellar age. Bontemps et al. (1996) made an important study of low-mass systems 
in Taurus and Ophiuchus, and found evidence for a secular decline in outflow power 
with age, while / remained constant. Intriguingly, the best correlation of outflow 
power was with circumstellar mass (as measured by the millimeter continuum flux) 
rather than with source bolometric luminosity; Saraceno et al (1996) also presented a 
similar correlation between millimeter continuum flux and outflow kinetic luminosity 
in systems with L < 10 3 L©. These results strongly suggest that the accretion rate 
and the outflow strength both decline in proportion to the disk and envelope mass. 

C. Deflected Infall Models 

Although the above analysis shows that / need not be necessarily higher in 
high-mass outflows, the very large masses involved, combined with the inefficiency 
of entrainment and momentum transfer by dense jets (especially once they escape 
their parent clouds), has led to some discussion of whether the sweeping up of 
ambient molecular gas by an accretion driven wind is a viable mechanism for these 
objects (e.g. Churchwell 1997a). 

A recent class of outflow models, termed circulation models, can naturally gen- 
erate outflow masses much greater than the stellar mass. In these, most of the 
infalling circumstellar material is diverted magnetically at large radii into a slow- 
moving outflow along the polar direction, while infall proceeds along the equatorial 
plane (Fiege and Henriksen 1996a, 1996b). The main attraction of these models is 
that they generate large outflow masses for even small stellar masses and can gen- 
erally explain the observed opening angles and velocity structure seen in high-mass 
systems. In particular, self-similar models predict that the velocity and density laws 
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should take the form V(r, 9) = U{6) ^GM^/r and p(r, 6) = /j(9)M*r- 3 (r/r ) 2a -°- 5 , 
where 0.25 > a > —0.5, r Q is an unspecified radial scale, and U(0) and ji{9) are 
dimensionless functions. It is then straightforward to show that the force and mass- 
flux in the outflow are related by F co /M C o = s/GM*/ R C o{U(0)} where (U(6)) 
= fU(0) 2 ii(6)duj/ $U(9)n(0)dw is the density-weighted average velocity over the 
outflow solid angle. The inferred (U(6)) (using the same M* as for our estimates 
of f v w /vk) is plotted in the top part of Fig. 5b as open circles. It is clear that 
high values (between 1000 and 10) are required. Generalized circulation models 
that include Poynting flux driving yield values of (U(6)) between 5 and 200 (Lery, 
Henriksen, Fiege, in preparation). In model cases where radiation transport is im- 
portant in setting up the flow, a power-law slope of 0.8 is predicted between Fco 
and Lboi, which is close to the observed slope ~ 0.7 (see also Henriksen 1994). Ob- 
servations seem to indicate a systematic decline of (U(8)) with Lboi, which would 
also have to be explained. 

There are several concerns about these circulation models. First, there are 
many unipolar CO outflows known, such as NGC2024-FIR5 (Richer et al. 1992), 
and the almost-unipolar HH46-47 system (Chernin and Masson 1991). These are 
naturally explained by swept-up wind models if the protostar is forming on the 
edge of a cloud or close to an HII region interface: the jet or wind propagating into 
the cloud will then sweep up a large CO flow, while in the opposite direction little 
evidence for a CO lobe will be seen. In circulation models, anisotropic solutions 
may also occur, but it is unclear why the weaker lobe would necessarily be on the 
side where the large-scale cloud density is low Second, in some objects such as B5 
(Velusamy and Langer 1998), there is an apparent lack of molecular material in the 
equatorial plane which can feed a circulation flow. Third, as discussed in section 
II, in some high-mass systems such as HH80-81 there is direct evidence for fast jets 
and bow-shock entrainment of molecular gas. Consequently, it seems more probable 
given the evidence presented that even high-mass outflows can be generated by the 
sweeping up of ambient cloud material by an accretion-driven stellar wind or jet. 
The details of the MHD driving mechanism in these cases, and of the momentum 
transfer between the wind and the jet, remain open issues. 



IV. SHOCK CHEMISTRY AND ENERGETICS 



The interaction between a supersonic protostellar wind and surrounding qui- 
escent material is expected to drive strong shock fronts. Shocks can be of type C 
(continuous) or J (jump) depending upon the shock velocity, the magnetic field, and 
the ionization fraction of the pre-shock gas (Draine and McKee 1993; Hollenbach 
1997). In the last few years, spectacular gains in sensitivity in the mm and IR do- 
mains have allowed us for the first time to witness the chemical and thermal effects 
of these shocks in molecular outflows. These observations yield direct estimates of 
the flow age, energetics, and entrainment conditions which represent an important 
new step toward a complete description of the outflow phenomenon. 
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A. Chemical Processing of ISM in Molecular Flows 

Theoretical expectations: By compressing and heating the gas, shock waves trig- 
ger new chemical processes which lead to a specific "shock chemistry" (see chapter 
by Langer et al., this volume). The most active molecular chemistry is expected to 
occur in C-shocks, as they increase the temperature to moderate values of about 
2000 K in a thick layer where molecules can survive, and reactions that overcome 
energy barriers can proceed. In particular, the very reactive OH radical can be 
formed by O+H2 — > OH + H (which has an energy barrier of 3160 K), and will 
contribute to the formation of H2O by further reaction with H2: OH + H2 — > H2O 
+ H (energy barrier: 1660 K). In dissociative J-shocks, molecules are destroyed in 
the hot (T ~ 10 5 K) thin post-shock layer and only reform over longer timescales, 
in a plateau of gas at ~ 400 K. Since some of these chemical processes are fast, 
and the cooling times are short, the chemical composition of the shocked regions is 
expected to be strongly time-dependent. 

Shocks also process dust grains. In the most violent J shocks, destruction of 
grain cores and thermal sputtering inject refractory elements (such as Si and Fe) 
into the gas phase (e.g. Flower et al. 1996). In slower C-type shocks, non-thermal 
sputtering will inject refractory and volatile species mainly from the grain mantles 
into the gas phase (e.g. Flower and Pineau des Forets 1994). The entrance of this 
fresh material, together with the high abundance of OH, will produce oxides such as 
SO and SiO (see Bachiller 1996, van Dishoeck and Blake 1998, and references there 
in). As the shocked gas cools, the dominant reactions will be again those of the 
usual low temperature chemistry, and depletion onto dust grain surfaces will reduce 
the abundances of some of the newly formed molecules (e.g. H2O, see Bergin et al. 
1998). However, the chemical composition of both the gas and the solid phases will 
remain altered with respect to pre-shock ones. 

Observations: The chemical effects of shocks have been observed in a number of 
outflows from low-mass Class objects, which are particularly energetic and contain 
shocked regions well separated spatially from the quiescent protostellar envelope. 
Recent examples include IRAS16293 (Blake et al. 1994; van Dishoeck et al. 1995), 
NGC1333 IRAS4 (Blake et al. 1995), and NGC1333 IRAS2 (Langer et al. 1996; 
Blake 1997; Bachiller et al. 1998). 

A comprehensive study of many different species has been recently carried out 
on L1157 (Bachiller and Perez Gutierrez 1997). The abundances of many molecules 
(e.g. CH3OH, H 2 CO, HCO+, NH 3 , HCN, HNC, CN, CS, SO, S0 2 ) are observed 
to be enhanced by factors ranging from a few to a few hundred. The extreme case 
is SiO which is enhanced by a factor of ~ 10 6 . There are significant differences 
in spatial distribution among the different species: some molecules such as HCO + 
and CN peak close to the central source, while SO and SO2 have a maximum in 
the more distant shocks, with OCS having the most distant peak. Other molecules 
such as SiO, CS, CH3OH, and H 2 CO show an intermediate behavior. Such differ- 
ences cannot be attributed solely to excitation conditions: an important gradient 
in chemical composition is observed along the outflow. It is very likely that this 
strong gradient is related to the time dependence of shock-chemistry. As an exam- 
ple, consider the chemistry of SO, SO2, and OCS, which has been recently modeled 
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by Charnley (1997). It is believed that sulfur is released from grains in the form of 
H2S, and that it is then oxidized to SO and SO2 in a few 10 3 yr. The formation of 
OCS needs a few 10 4 yr. This is in general agreement with observations, since the 
SO/H2S and SO2/H2S ratios do increase with distance from the source (i.e. with 
time), and OCS emission is only observed in the most distant position (i.e the oldest 
shock). Hence, chemical studies are of high potential to constrain the age and time 
evolution of molecular outflows. 

B. Shock Cooling and Energetics 

Millimeter observations of shock-enhanced molecules trace chemically processed 
gas that has already cooled down to 60-100 K, as indicated e.g. by multi-line 
NH 3 studies (Bachiller et al. 1993, Tafalla and Bachiller 1995). Emission from 
hotter post-shock gas, on the other hand, is important to obtain information on the 
instantaneous energy input rate and pre-shock conditions in outflows. 

Hot (T > 1000 K) shocked gas: Hollenbach (1985) suggested that the [O I] 63/xm 
line should offer a useful, extinction-insensitive measure of dissociative J-shocks in 
molecular flows. Because [O I] 63/um is the main coolant below ~ 5000 K, its 
intensity is roughly proportional to the mass flux into the J-shock, Mjs, through 
the relation L^qj^/Lq = 10 4 x Mjs/Mqjt -1 , as long as the line remains optically 
thin (i.e. n a Vjs < 10 7 km s _1 cm -3 , where n a is the pre-shock density and Vjs is 
the J-shock velocity). 

First detections of [O I] 63/um in outflows were obtained with the Kuiper Air- 
borne Observatory (KAO) toward 3 HH objects and 5 highly collimated Class 
outflows (Cohen et al. 1988; Ceccarelli et al. 1997). Since the advent of the In- 
frared Space Observatory (ISO), the Long Wavelength Spectrometer (LWS) has 
revealed [O I] 63/um emission in at least 10 more HH objects and molecular out- 
flows (Liseau et al. 1997; Saraceno et al. 1998). These authors find a surprisingly 
good correlation between current values of Mjs derived from [O I] 63/um and time- 
averaged M a ve values derived from mm observations of the outflow assuming ram 
pressure equilibrium at the shock (i.e. Fco = M ave x Vjs)- Both values agree for 
Vjs ~ 100 km s _1 . The dispersion in this correlation, roughly a factor of 3, is of 
the same order as the uncertainties in Fqo caused by opacity and projection effects 
(Cabrit and Bertout 1992). Hence, CO-derived momentum rates in outflows do not 
appear to suffer from large systematic errors. 

With the development of large format near-IR arrays in the early 1990's, it has 
also become possible to map molecular outflows in the 2.12 / um v = 1 — S(l) line of 
H2, a tracer of hot (~ 2000 K) shocked molecular gas. In both low-luminosity and 
high-luminosity outflows, the H2 emission delineates single or multiple bow-shaped 
features associated with the leading edge of the CO emission (Davis and Eisloffel 
1995; Davis et al. 1998), as illustrated in Fig. 1 for HH211; in a few cases, H2 
emission also traces collimated jets and cavity walls (e.g Bally et al. 1993; Eisloffel 
et al. 1994). Observed surface brightnesses and rotational temperatures ~ 1500- 
2500 K indicate moderate velocity shocks: either J-shocks of speed 10-25 km s _1 or 
C-shocks with V s ~ 30 km s _1 and low filling- factor (Smith 1994; Gredel 1994). In 
particular, the morphology, line profile shapes, intensity, and proper motions of H2 
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2.12^111 bows are well reproduced by hydro dynamical simulations of jets propagating 
into the surrounding cloud, where H2 2A2fxm emission arises mostly in the non- 
dissociative wings of the bowshock (Raga et al. 1995; Micono et al. 1998; Suttner 
et al. 1998). 

If these bowshocks are also where most of the slow molecular outflow is being 
accelerated, and if H2 emission dominates the cooling (as expected e.g. in 2000K 
molecular gas at densities of 10 5 — 10 8 cm -3 ), then L(H2)/Lco should be of order 
unity (see e.g. Hollenbach 1997). In the 5 flows studied by Davis and Eisloffel (1995), 
the observed ratio L(H2)/Lco has a median value of ~0.4, but it covers a very broad 
range from 0.001 to 30. The discrepancies could be caused by uncertainties in 2/im 
extinction (corrections typically amount to 10-100 for Ay = 20 — 50 mag), by the 
use of unreliable Leo estimates, or — in the case of very low ratios — by a strong 
decrease in outflow power over time (see W75N; Davis et al. 1998). Hence the H2 
2.12//m line alone is not a sufficient diagnostic of the outflow entrainment process. 

Warm (T ~ 300-1000 K) molecular gas: The ISO mission has led to the de- 
tection of a new component of warm post-shock gas at T ~ 300-1000 K in several 
outflows. Fig. 6 shows a map of the L1157 outflow in the v = — S(5) pure rota- 
tional line of H2 at 6.9/um obtained with ISOCAM (Cabrit et al. 1998). A series of 
bright emission spots are seen along the outflow axis. They coincide spatially with 
hot shocked gas emitting in the H 2 2.12 / um line (Eisloffel and Davis 1995), and with 
the various peaks of shock-enhanced molecules identified by Bachiller and Perez 
Gutierrez (1997; see Sect. IV. A). However, they trace an intermediate temperature 
regime of ~ 800 K, considerably lower than the 2000 K observed in ro-vibrational 
H 2 lines. Warm H 2 at 700-800K was also found with the ISO Short Wavelength 
Spectrometer (SWS) in two outflows from very luminous sources, Cep A and DR 21 
(Wright et al. 1996; Smith et al. 1998). Finally, warm CO at T ~ 330 - 1600LT 
was detected in high- J lines ( J up = 14 to 28) with ISO-LWS toward 5 outflows of 
various luminosities, while H2O and OH lines were detected in two cases (Nisini et 
al. 1996, 1997; Ceccarelli et al. 1998). 

The observed emission fluxes and temperatures in H 2 and CO are well explained 
by non-dissociative J-shocks with V s ~ 10 km s -1 or slow C-shocks with V s ~ 10- 
25 km s -1 , and n a ~ 10 4 — 3 x 10 5 cm -3 (e.g. Wright et al. 1996; Nisini et al. 
1997; Cabrit et al. 1998). One important constraint is the rather low [H 2 0]/[H 2 ] 
abundance ratio ~ l-2x 10~ 5 observed in HH54 and IRAS16293 (Liseau et al. 1996; 
Ceccarelli et al. 1998); steady-state C-shocks would predict complete conversion of 
O into water. The relatively high [OH]/[H20] ratio ~ 1/4-1/10 in these two flows 
suggests that the shock age is too short for conversion to be complete, and points to 
the need for time-dependent C-shock models for proper interpretation of the data 
(e.g. Chieze et al. 1998). 

Mid and far-infrared emission from this warm molecular gas component appears 
more tightly correlated with Leo than the 2/im H2 lines. In 4 out of the 5 outflows 
studied by Nisini et al. (1997), the FIR CO luminosity represents 10-30% of the 
flow kinetic luminosity, in good agreement with C-shock calculations in the inferred 
density and velocity range (Kaufman and Neufeld 1996). The only large discrepancy 
is observed in IC1396N, an object contaminated by PDR emission (Molinari et al. 
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1997). In L1157, the H 2 luminosity of warm gas is also around 10% of Leo (Cabrit 
et al., 1999). Hence, these slow shocks seem sufficient to drive the whole outflows. 
Detailed comparisons between H2 and FIR-CO lines in the same objects are now 
under way to further narrow the range of possible shock models and perhaps allow 
us to discriminate between wide-angle wind and jet scenarios for the entrainment 
of outflows. 
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Figure 6. The LI 157 outflow mapped in the 6.9/im pure rotational line of H2 
(adapted from Cabrit et al. 1998) with CO (2-1) contours superimposed (from Bachiller 
and Perez Gutierrez 1997). 
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V CONCLUSIONS 

We have shown that the current data on the structure and energetics of molec- 
ular outflows suggest broad similarities across the entire luminosity range, from 1 to 
10 5 L©. If the flows are swept up by a stellar wind or jet, we find that M w v w / M a VK 
has a value of about 0.3 for all flows, perhaps suggesting that flows have a common 
drive mechanism. However, it remains unclear if the MHD disk and X-wind mod- 
els which have been used to explain low-mass outflows are appropriate in the very 
different physical regime of high-mass YSOs. While observational data continue 
to improve these constraints, there remains an urgent need for a larger sample of 
molecular outflows, particularly from high-mass stars, to be fully mapped at high 
resolution; at the moment it is very possible that our estimates of the properties 
of high-mass systems are biased by strong selection effects. Single dish data, es- 
pecially from the new focal plane arrays, plus interferometric images at millimeter 
wavelengths will continue to accumulate. However, only when the large millime- 
ter interferometer (MMA/LSA) is operational will it be possible to acquire high- 
resolution data quickly enough to study large samples of outflows in detail. 

The nature of the shocks which drive outflows is slowly becoming clearer. We 
now have diagnostics of all the temperature components in the outflows, from the 
2000K gas seen in the 2 fim H2 lines, to the several hundred Kelvin component 
recently detected by ISO, through to the cool massive component seen in the mil- 
limeter waveband where most of the momentum is eventually deposited. The rela- 
tionship between these components is providing valuable tests of the outflow mech- 
anism, although a fuller understanding will require observations at higher angular 
resolution than ISO provided. The SOFIA and especially the FIRST missions will 
provide valuable data in this area. 
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